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An effective meson Lagrangian including a scalar glueball is constructed 
on the base of U{2>) x f7(3) chiral symmetry. The glueball is introduced into 
the meson Lagrangian by using the principle of scale invariance of an effective 
Lagrangian and the dilaton model. The singlet-octet mixing of scalar meson 
states is described by means of 't Hooft interaction. The contribution of 
the scalar and pseudoscalar anomalies into the breaking of scale invariance is 
taken into account. The mixing of quarkonia with the glueball is described. 
The mass spectrum of scalar mesons together with the glueball and also their 
strong decay widths are calculated. Prom comparing the obtained results 
with experimental data, it follows that /o(1500) is rather a glueball, whereas 
/o(1710) is a quarkonium. This accords with the results obtained in our 
previous works where radially-excited scalar meson states were described. It 
is shown that p-mesons play an important role in the description of glueball 
decays. 
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1. INTRODUCTION 



The self- interaction of gluons, a peculiarity of QCD, gave an idea that 
gluons can form bound states that can propagate as particles in the space. 
Unfortunately, because of theoretical problems, there is no yet the exact an- 
swer to whether these states really exist or not. However, from recent lattice 
simulations one can conclude that it is most probably that glueballs are 
real objects of our world. Having assumed that glueballs exist, one can try 
to construct a model to describe their interaction with other mesons, their 
properties, such as, e. g., mass and decay width, and to identify them with 
observed resonances. 

An exact microscopic description of bound gluon states cannot be done 
systematically in the framework of QCD. In this situation, QCD-motivated 
phenomenological models are the tool that can help to deal with glueballs 
as well as with quarkonia which form the most of observed meson states. 
However, using these models to describe glueballs, we encounter many dif- 
ficulties concerning, e. g., the ambiguity of the ways of including glueballs 
into models and identification of experimentally observed meson states. This 
explains the variety of points of view on this problem. 

First of all, we do not know the exact mass of a glueball. From the 
quenched QCD lattice simulations, Weingarten (see, e. g., |T],0]) concluded 
that the lightest scalar glueball is expected around 1.7 GeV. Amsler [Q] con- 
sidered the state /o(1500) as a candidate for the scalar glueball. QCD sum 
rules []5j] and K-matiix method showed that both /o(1500) and /o(1710) 
are mixed states with large admixture of the glueball component. 

All bound isoscalar qq states are subject to mixing with glueballs, and 
their spectrum has many interpretations made by different authors. For in- 
stance, Palano [|7| suggested a scenario, in which the states ao(980), Kq{14:30), 
/o(980), and /o(1400) form a nonet. The state /o(1500) is considered as 
the scalar glueball. Tornqvist et al. [jSj looked upon the states /o(980) and 
/o(1370) as manifestations of the ground and excited ss states, and upon 
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the state /o(400 — 1200) as the ground uu state. Eef van Beveren et al. 
considered the states /o(400 — 1200) and /o(1370) as ground uu states, and 
the states /o(980) and /o(1500) as ground ss states. Two states for each 
qq system occur due to pole doubhng, which takes place for scalar mesons 
in their model. Shakin et al. [ pXI| , pr[|] obtained from a nonlocal confinement 
model that the /o(980) resonance is the ground uu state, and /o(1370) is the 
ground ss state. The state /o(1500) is considered as a radial excitation of 
/o(980). They believe the mass of scalar glueball to be 1770 MeV. 

In our recent papers []T2[], following methods given in [[T3|-[ro|, we showed 



that all experimentally observed scalar meson states with masses in the inter- 
val from 0.4 to 1.71 GeV can be interpreted as members of two scalar meson 
nonets — the ground state of the meson nonet (lighter than 1 GeV) and its 
first radial excitation (heavier than 1 GeV). We considered all scalar mesons 
as qq bound states and took into account the singlet-octet mixing caused by 
't Hooft interaction. In [|T2|, we obtained a scalar isoscalar state with mass 



1600 MeV and had to choose to which of the experimentally observed states, 
/o(1500) or /o(1710), we should ascribe it. From our analysis of the strong 
decay rates calculated in our model we found that /o(1710) fits to the nonet 
of quarkonia better than /o(1500). Therefore, we supposed that the state 
/o(1500) contains greater admixture of the scalar glueball (see [@,|5j). How- 
ever, the final decision should be made after including the scalar glueball into 
the model and taking account of its mixing with quarkonia. In the present 
work, that is devoted to solving this problem, from the analysis of strong 
decay widths of the glueball we again come to an analogous conclusion]^. 

To describe the properties of the glueball and its interaction with quarko- 
nia, one should introduce a scalar isoscalar dilaton field x i^ito our model, 
in addition to the quarkonia that have already been described []T^] . For this 



purpose, one can make use of the idea of approximate scale invariance of 



However, radially-excited states are not yet considered. 
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effective Lagrangians based on tfie dilaton model. Such models were studied 
by many autfiors (see, e. g., [[r7|-PT|])- Unfortunately, there is no unique way 
to introduce the dilaton field into a chiral Lagrangian. This explains the 
large number of models dealing with glueballs. 

The guideline one should follow when introducing the dilaton field into 
an effective meson Lagrangian is to reproduce the Ward identity connected 
with the scale anomaly. The latter leads to the following equation for the 
vacuum expectation value of the divergence of the dilatation current: 

{d,S^)=Cg- E rn'^im). (1) 

q=u,d,s 

where is the number of colours; Nf is the number of flavours; C^j^) and 
{qq) are the gluon and quark condensates; is the current quark mass. 

In this paper we are going to use the most natural method of introduc- 
ing the dilaton field into the effective Lagrangian by requiring that, in the 
chiral limit, our Lagrangian should be scale-invariant except for the dilaton 
potential and terms induced by gluon anomalies. To realize this program, 
one should multiply all dimensional parameters of the original Lagrangian 
(without dilaton) by a corresponding power of the dilaton field divided by 
its vacuum expectation value Xc- Thus, instead of the four-quark coupling 
constant G, the 't Hooft couphng constant K, ultraviolet cutoff A (necessary 
for regularizing the divergent integrals coming from quark loops), and the 
constituent quark masses niq {q = w, s), one should use G'(Xc/x)^? ^(Xc/x)^? 
Mx/XcX and mq{x/Xc)- 

Current quark masses are not multiplied by the dilaton field and violate 
scale invariance explicitly, as it takes place in QCD. Their contribution to 
the divergence of dilatation current is determined by quark condensates and 
disappears in the chiral limit (see (|T])). 

The scale invariance is also broken by those terms in the effective La- 
grangian that are induced by the pseudoscalar and scalar gluon anomalies 
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and look as follows [|22|,|23 



^an = -/i00O + (3) 



where h^.hcj are constants. 0o = \/2/30if — a/1/30s, 0o and do ((c^o) 7^ 0) 



do = — ^XjZos^ where 7^ 0) consists of w((i)-quarks and Og 

{{(Js) 7^ 0) of s-quarks. 

These terms appear due to the 't Hooft interaction. When restoring scale 
invariance of the effective Lagrangian by inserting dilaton fields (the proce- 
dure of the restoration of scale invariance is given in Sect. 3), these terms 
must be treated separately. Moreover, it turns out that these terms deter- 
mine the most of quarkonia-glueball mixing. 

Omitting, for a moment, the 't Hooft interaction in our approach, we 
require the Lagrangian to be scale-invariant in the chiral limit both before 
and after the spontaneous breaking of chiral symmetry (SBCS), except for 
the dilaton potential. This property can be obtained by considering (after 
bosonization when the effective Lagrangian is expressed in terms of bosonic 
scalar and pseudoscalar fields a and 0) the shift of the scalar meson field a 

a = a'-m^, (m^ = 0), (4) 

Xc 

where (o"')o = 0, {a)o = — m, guaranteeing that the relation (|I]) is satisfied. 
The nonzero vacuum expectation value of a appears as a result of SBCS, and 
thus, the constituent quark mass m is produced. In the case of nonvanishing 
current quark masses, (g) changes by including an additional (nonscaled) 
mass term mP into the r.h.s. 

a = a — m— + m^. (5) 

Xc 

The structure of the paper is as follows. In Section 2, we derive the usual 
U{3)xU (3)-flavour symmetric effective Lagrangian with the 't Hooft interac- 
tion and without dilaton fields. In Section 3, the dilaton field is introduced 
into the effective Lagrangian obtained in Section 2. In Section 4, the gap 
equations are investigated, the quadratic (in fields) terms are deduced and 



the mixing matrix for scalar isoscalar states is introduced. In Section 5, the 
numerical estimates for the model parameters are given. The main strong 
decays of scalar isoscalar mesons are calculated in Section 6. It is shown 
there that p-mesons play an important role in the decay of a glueball into 
four pions. Finally, in the Conclusion, we discuss the obtained results. 

2. CHIRAL EFFECTIVE LAGRANGIAN WITH 'T HOOFT INTERACTION 

A U{3) X U{3) chiral Lagrangian with the 't Hooft interaction was inves- 
tigated in paper It consists of three terms (see below). The first term 



represents the free quark Lagrangian, the second is composed of four-quark 
vertices as in the NJL model, and the last one describes the six-quark 't 
Hooft interaction []25[] that is necessary to solve the Ua{^) problem. 



L = q{zd - m% + f E [{qKqf + {qil.Kqf] - 

-K {det[g-(l + 75)g] + det[g-(l - 75)^]} . (6) 

Here G and K are coupling constants, Aa (a = 1,...,8) are the Gell-Mann 



matrices, Aq = y^2/3 1, with 1 being the unit matrix; rnP is a current quark 
mass matrix with diagonal elements m^, m[), (m^ « m[|). 

The standard bosonization procedure for local quark models consists in 
replacing the four-quark vertices by Yukawa coupling of quarks with bosonic 
fields which enables one to perform the integration over quark fields. The 
final effective bosonic Lagrangian appears then as a result of the calculation 
of the quark determinant. To realize this program, it is necessary, using the 
method described in [j2^-|27|], to go from Lagrangian (p) to an intermediate 



Lagrangian which contains only four-quark vertices 

n 1 ^ (—'] (-1.'] 

L = q{id -m)q+- [Gab (Q'^aq) {qnq) + G\f, {qilhTaq) [qilbUq)] , (7) 

^ a,h=\ 

where 
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ra = Xa (a = l,...,7), rg = (v^Ao + A8)/V3, 

rg = (-Ao + V2A8)/V3, 

G^n^ = = Gg^ = G± AKmjtim,), 

Git^ = = = = G ± 4irm,/i^(m0, 

Gg) = G T 4A^m,/i^(m,), G^^ = G, 

Cst^ = G'g^ = ±4v^A^m,/f (m,), 

Git^ = {a^h- a,6=l,...,7), 

^is'' = = Gg"* = Gga'' = (a = 1, . . . , 7), (8) 

and is a diagonal matrix composed of modified current quark masses: 

ml = ml- ?>2KmumJ^{mu)I^{m,), (9) 
jnl = ml- ?>2Kmll^{muf. (10) 

Here niu and are constituent quark masses and the integrals 

In{ma) = ^J dlk ^}f~^2 ^ (n = l,2; a = n,s), (11) 



are calculated in the Euclidean metric and regularized by a simple 0(4)- 
symmetric ultraviolet cutoff A. For /^(ma) one gets 

AVj^^(A^-™^ln(A! + l]], (12) 

where nia represents a corresponding constituent quark mass: or mg. Note 
that we have introduced the notation of constituent quark mass already here, 
although they will be consistently considered only later, when discussing mass 
gap equations (compare ( |5tl|) and ([5T|) below) and the related shift of scalar 
meson fields. However, as we want to use an effective four-fermion interaction 
instead of the original six-quark one, we have to use full quark propagators 
with constituent quark masses to calculate quark loop corrections for the 
constant G (see (g)). For the definition of the constituent quark masses see 
(Ol) and (^ below. 
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In addition to the one-loop corrections to the constant G at four-quark 
vertices, we modified the current quark masses m[| (see (H) and (|TDD). This 
is to avoid the problem of double counting of the 't Hoot contribution in gap 
equations which was encountered by the author in [j27[] . After the redefinition 
of the constant G and of the current quark masses, we can guarantee that in 
the large- A^c limit the mass spectrum of mesons and the gap equations, de- 
rived from the new Lagrangian with modified four-quark vertices and current 
quark masses, are the same as those obtained from the original Lagrangian 
with six-quark vertices. 

Now we can bosonize Lagrangian ([7]). By introducing auxiliary scalar a 
and pseudoscalar fields, we obtain [[r3| , [n| , pi|] 



-^Trln|z^-mO+ |]r„K + z75(/)a)|. (13) 

As we expect, the chiral symmetry is spontaneously broken due to strong at- 
traction of quarks in the scalar channel and the scalar isoscalar fields acquire 
nonzero vacuum expectation values (cra)o 7^ (a = 8,9). These values are 
related to basic model parameters G, m^, and A via gap equations as it will 
be shown in the next Section. Therefore, we first have to shift the a fields by 
a proper value so that the new fields have zero vacuum expectation values: 

(ya = (y'a- ^a^ A (^a)o = 0, (14) 

where /ia = 0, (a = 1, . . . , 7), /ig = rriu, /ig = —ms/V^ and ji^ = 0, (a = 
1, . . . , 7), /ig = m^, /ig = — After this shift we obtain: 

/:(cr',0)= LG(cr', 0) Tr In I ?a-m+ E '^a(c^l + ^75 0a) I , (15) 



where 



1 ^ -1 

^ a,6=l 

-I E 0a(G'W);>,, (16) 
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and m is a diagonal matrix of constituent quark masses for different flavors. 
From Lagrangian (|T5D we take only those terms (in momentum space) which 
are linear, squared, cubic, and quadruple in scalar and pseudoscalar fields-f^ 

C{u', 0) = LG(a', 0) + tr I^{m){{d^,a'f + {d^c\>f) - ^ml^{m)a' + 
+ 2/f (m)(a'2 + 02) _ 4m2/2^(m)a'2 + 



+ 4m/2^(m)a'(a'2 + <\?f - I^{m){a'^ + 02)2 + 



+ /2^(m)[a'-m,0]: 

9 

a=l 



= J2 (f^aTa: 
a=l 



(17) 
(18) 



where "tr" means calculating the trace over r-matrix expressions and [...]- 
stands for a commutator. (The calculation of "tr" is explained in details in 
T3|] ) The expression for /{^(ttt-q) in Euclidean metric is given in ([T^). The 



integrals /2^(ma) are also calculated in Euchdean space-time 



In I At + 1 



A2 



A2 + m 



(19) 



Then, we renormalize the fields in ( \n\) so that the kinetic terms of the effec- 
tive Lagrangian are of the conventional form, and diagonalize the isoscalar 
sector. 



"1 

4 



+tr 



+\[m, 0^- - m^<j' ' + mga'{a' ^ + Z0^ ^) - ^[m, 01- 
-^{{a'^ + Z(f)'^f-[a\(J)']l) 



r \ ^ r 



(20) 
(21) 



Despite that the scalar fields are of the main interest in this paper, we still need pseudoscalar fields to 
fix the model parameters. 
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For Lq we have: 



2 a,6=l 



a6 



Here we introduced Yukawa coupling constants 



(J. 



= '/Zga(j)'a, 



(22) 



(23) 



2 2 2 2 2 F/i rA/ 
fi'l = 5^2 = 5^3 = 5'8 = = [4^2 ("^ 

fi'4 = 5^5 = 5^6 = 5'7 = [4/^("^w, "^s)]~\ 

^(A2 - A;2) 



2{mu,ms] 



_Nc 
(27r) 



d k— 



(47r)2(m^ — m 



2. / A' 
m„ In I — - + 1 



m^, In 



mi 



Z= 1 



Gm^ 



1.44. 



(24) 



A2 



+ 1 



(25) 
(26) 



where we have taken into account 7r-74i-transitions leading to an additional 



Z factor, with Ma^ being the mass of axial- vector meson (see []T3|). The 



renormalized scalar and pseudoscalar fields in ([2D|) -(P3|) are marked with the 
superscript r. 

The mass formulae for isovectors and isodublets follow immediately from 
( pO|) . One just has to look up the coefficients at a^^ and 0^^. There are 
still nondiagonal terms in ( P2| ) in the isoscalar sector. This problem is solved 
by choosing the proper mixing angles both for the scalars and pseudoscalars 
(see, e. g., [)2^). As we are going to introduce the glueball field, the mixing 
with scalar isoscalar quarkonia will change the situation. One has to consider 
the mixing among three states, which cannot be described by a single angle. 
For simplicity, in our estimations we resort to a numerical diagonalization 
procedure, not to the algebraic one. Concerning the pseudoscalar sector, one 
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can avail oneself with the results given in []2^. All what concerns dealing 
with the glueball is discussed in the next Section. 



3. NAMBU-JONA-LASINIO MODEL WITH DILATON 

As we have already mentioned above, we introduce the glueball field into 
our effective Lagrangian obtained in the previous Section, as a dilaton. For 
this purpose, we use the following principle. Insofar as the QCD Lagrangian 
is scale-invariant in the chiral limit, we suppose that our effective meson 
Lagrangian, motivated by QCD, has also to be scale- invariant both before 
and after SBCS in the case when the current quark masses are equal to 
zero. As a result, we come to the following prescription: the dimensional 
model parameters G, A, and rUa are replaced by the following rule: G 

G{Xc/x)^, K K{xc/xf, A A(x/Xc)^ rUa ma{x/Xc), where x is the 
dilaton field with the vacuum expectation value Xc- But there are terms that 
break scale invariance. They are the terms containing current quark masses; 
the scale anomaly of QCD, reproduced by the dilaton potential; and terms of 
the type /i^^o ^^^1 hadQ (see (|3D), induced by gluon anomalies in the meson 
Lagrangian. 

As it was mentioned in the previous paragraph, the current quark masses 
break scale invariance and, therefore, should not be multiplied by dilaton 
fields. The modified current quark masses fnP are also not multiplied by 
dilaton fields. In particular, this transforms formula (|T^) to what follows: 



aa = a'^- Ha— + f^t (27) 

Xc 

Finally, we come to the following Lagrangian: 

£(a^ 0^ x) = c{x) + i^kin(a^ 00 + ^G(a^ 0^ x) + ^^i-ioop + al^h. (28) 

Here C{x) is the pure dilaton Lagrangian 

^X) = lid.xf - VU) (29) 
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with the potential 



V{x) = B (^) 
VXo/ 



In ^ 



1 



(30) 



that has a minimum at x = Xo? and the parameter B representing the vacuum 
energy when there are no quarks. The curvature of the potential at its 
minimum determines the bare glueball mass 

A^^B 



nig = 



xo 



(31) 



The part Lkin(cr'', (jf) of Lagrangian (|2^) contains pure kinetic terms: 

^kin(a^ 00 = ^ E [{d,<f + {d,¥af) ■ (32) 
The next term reads 



^g(^O0Ox) 

2 



\XcJ nb^l V Xc / ' 



-1 

ab 



X 



X ( ^5(7^ - h /i^ 

V Xc 
2 9 



^ VXc/ o,5=l 



(33) 



The sum of one-loop quark diagrams is denoted as Li_ioop 

^3 

-^i-ioop = tr 

2 



VXc/ 

22r2/X\ I X r/'r2,/r2\ 

VXc/ Xc 



.Xc 
2\2 



(34) 



As one can see, expanding ( |33|) in a power series of we can extract 
a term that is of order x^- It can be absorbed by the term in the pure 
dilaton potential which has the same degree of x- Obviously, this leads only 
to a redefinition of constants B and xo which anyway are not known from 
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the very beginning. Moreover, saying in advance, the terms hke do not 
contribute to the divergence of the dilatation current (|I]) because of their 
scale invariance. 

If the procedure of the scale invariance restoration of this Lagrangian 
is implemented, the part induced by gluon anomalies also becomes scale- 
invariant. To avoid this, one should subtract this part in the scale-invariant 
form and add it in a scale-breaking (SB) form. This is achieved by including 
the term ALan: 

ALan=-Lanf-)' + Lf. (35) 
\XcJ 

The term Lan was introduced in (|^). In Lan, we will use renormalized ctq 
and (f)'^ fields in place of ctq and 00, however, taking into account the effects 
of nonzero vacuum expectation value of ctq. Let us define the scale-breaking 
term L^^. The coefficients hg- and /i^ in (0) can be determined by comparing 
them with the terms in ([531) that describe the singlet-octet mixing. We obtain 



3 3 

K = -^^9u9sZ{G'^^^)^i, K = :^9ugs{G^~\^- (36) 

If these terms were to be made scale invariant, one should insert {x/Xcj^ i^ito 
them. However, as the gluon anomalies break scale invariance, we introduce 
the dilaton field into these terms in a more complicated way. The inverse 
matrix elements (G'*^^-')^^^ and {G'^^^^^^ 



ab 1 

^88 ^99 ~ l^g" 



^89 - (-) (-) .^(-)^2' ^-^^^ 
^88 ^-^99 ~ I'J'89 ; 

are determined by two different interactions. The numerators are fully de- 
fined by the 't Hooft interaction that leads to anomalous terms (0) breaking 
scale invariance, therefore we do not introduce here dilaton fields. The de- 
nominators are determined by constant G describing the main four-quark 
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interaction, and the dilaton field is inserted into it, according to the pre- 
scription given in the beginning of this Section. Finally, we come to the 
following form of L^^: 

= (- h,4>I> ' + K (c^S -F,>^^ + F») ' j {£j\ (39) 
From it, we immediately obtain the term ALan: 

ALan= (MS'-^.f^S-^0- + Fo"V| (-1 I- (41) 



Fo = + F« = + (40) 





Let us now consider the vacuum expectation value of the divergence of the 
dilatation current calculated from the potential of the effective meson-dilaton 
Lagrangian: 

( ^ dV dV 



X=Xc 



4 



= Ab{^\ -2K{F,-F^f - E Ki^q)- (42) 



Here V = V{x) + ^(c^^? 0^^? x)? S'lid V{(t^ , (jf , x) is the potential part of La- 
grangian C{(T^ ^ (jf ^ x) that does not contain the pure dilaton potential. The 
expression given in (|42| ) is simplified by using the following relation of the 
quark condensates to integrals Ii{mu) and /^(mg): 



Amql^iniq) = -{qq)o, {q = u,d,s), (43) 
and that these integrals are connected with constants ^ through gap equa- 



tions, as it will be shown in the next Section (see (^) and (^) below). Com- 
paring the QCD expression (|I]) with (|2|), one can see that the term E rn^qiqq) 
on the left-hand side is canceled by the corresponding contribution on the 
right-hand side. Equating the right hand sides of (|I]) and (^), 
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Cg- E mq{qq) = 



q—u,d,s 



^XO/ q~u,d,s 



we obtain the correspondence 

4 9 



(44) 



Cg = 45 m + E (a^S - f^'a){G^-YJ{fi, - fit) - 2K {Fo - n 

VXO/ a,6=8 



0\2 



(45) 



where = (a = 1, . . . 7), ^Ug = m^, and /ig = —ms/V^. This equation 
relates the gluon condensate, whose value we take from other models (see, 
e. g., [^), to the model parameter B. The next step is to investigate the 
gap equations. 



4. EQUATIONS 

As usual, gap equations follow from the requirement that the terms linear 
in and x' should be absent in our Lagrangian: 



5C 



6al 



= 0, 



6jC 



= 0. 



(46) 



Here, the field x' = X ~ Xc with a zero vacuum expectation value {x')o = 0, 
is associated with the glueball field. In further calculations, the Lagrangian 
is expanded in power series of x'- As a result the following equations are 
obtained: 

(m. - rnl){G^-\i - "^^^^{G^-^^ - ^mj^{m^) = 0, 



V2 

{rris - fhl){G^-\^ - V2{mu - fril){G^-\l - ^m,I^{m,) = 0, 
AB (^)' - In (^) V 1 f E nl{G^-\M - 3/i5)) - 

VXo/ XO VXo/ Xc \a,b^8 J 



(47) 
(48) 



2h 



Xc 



- {Fo - F, 



0\2 




0. 



(49) 



Using (0) and ([TDD, can rewrite equations ( |47|) and (|4^) in the well- 
known form [ P7|] : 
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= m„ - 8Gm„/f(m„) - 32Km^mJi(m^)I^(m,), 
m° = m, - 8Gm.Ii(m,} - 32K(mJi(m^}f. 



(50) 
(51) 



To define the masses of quarkonia and the glueball, let us consider the 
part of Lagrangian (^) which is quadratic in fields and x' which we 
denote as L^^^ 



9899 



— )\ — 1 r r 
/89 ^8^9 



+ [G )ab9b(^bX + - 



a, 6=8 Xc 



Xc 



V3 



x'. 



(52) 



where 



+ E 4/iO(G^-0„t(/i6-/ii:)-MF(? + 4/..Wn 

a,b=8 



a,b=i 



i0\2> 



(53) 



is the glueball mass before taking account of mixing effects. 

From this Lagrangian, after diagonalization, we obtain the masses of three 
scalar meson states: ai, an, and ani, and a matrix of mixing coefficients b 
that connects the nondiagonalized fields = cr^, erg = ds, x' with the physical 
ones di, dii, diii 



fill 



(54) 



V ^x'fTi ^x''^!! K'^m J \ ^ni / 



5. MODEL PARAMETERS AND NUMERICAL ESTIMATES 



The basic parameters of our model are G, A, m^, and m^. After the 
dilaton fields are introduced, they keep their values [P^ , [2^] : 
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m 



280 MeV, m. = 420 MeV, A = 1.26 GeV, 



G = 4.38 GeV~2, K = 11.2 GeV^^ (55) 

Moreover, new three parameters xo? Xc? cind B appear. To fix the new 
parameters, one should use equations (|451) , (1491) , and the physical glueball 



mass. As a result we obtain for xo a-^d B: 

'""^ ' 4[C, - (/i2 - /i2)(G(-)),,i(/i5 - nl) + 2/.. (Fo - F^,f] ' ' ^ ^ 

= i • ^''^ 

We adjust the parameter Xc so that the mass of the heaviest scalar meson, 
(Tin, would be either 1500 MeV or 1710 MeV. The result of our fit for both 
cases is given in Table One also will find the mixing coefficients in Table 0. 



6. DECAY WIDTHS 



Once all parameters are fixed, we can estimate the decay widths for the 
main strong decay modes of scalar mesons: ai t^tt^KK, ijij, ijij', and 47r 
where / = I, II, III. 

Note that, in the energy region under consideration (~ 1500 MeV), we 
work on the brim of the validity of exploiting the chiral symmetry that was 
used to construct our effective Lagrangian. Thus, we can consider our results 
as rather qualitative. 

Let us start with the lightest scalar isoscalar meson state di, associated 
with /o(400 — 1200). This state decays into pions. This is the only strong 
decay mode, because ai is too light for other channels to be open. The 
amplitude describing its decay into pions has the form: 

^.1^.+.- = '2A'^^^-b^'a, + 2^"6,„,,, (58) 
Ai,,- = A- = 2g,m,Z, (59) 

Xc 
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where A'^^+^- is the contribution from the glueball component; and ^4", from 
the {uu) quarkonium one. The coefficients 6^^! and b^^ai represent the corre- 
sponding elements of the 3x3 mixing matrix for scalar isoscalar states (see 
Table g). Both contributions have equal signs and add to the the width of 

To calculate the decay width of a meson into two mesons, one can use the 
following formula: 

\A\' \'l\M\MlMl) 
^ = TfeM3 . (60) 

where A is the amphtude of the process; M is the mass of a decaying particle; 
Ml and M2 are masses of secondary particles; r is the dimension of the 
permutation symmetry group in the phase space of final states. The function 
X{x,y,z) is defined as follows [pD|] : 



X{x,y,z) = {x - y - zf - 4yz. (61) 

For the decay of ai into pions, formula ( |60|) can be rewritten in a simpler 
form 



1 - (62) 



Using isotopic symmetry, we obtain the total width 

r^i^.. = ^r,j^,+,- « 820 MeV (63) 

for the case when the model parameters are fixed for the state din identified 
with /o(1500), and 

r^j^^^ « 830 MeV. (64) 

for the case am = /o(1710). The experimental value is known with a large 
uncertainty and is reported to lie in the interval from 600 to 1000 MeV ^9 



The amplitude describing the decay of the state an which we identify with 
/o(980) into pions also consists of two parts 



Here the glueball contribution is small again and the quarkonium determines 
the decay width, however, in this case both contributions are opposite in sign 
and slightly compensate each other. The width of the state an is close to 
that obtained in the model without glueballs . We obtain 



if (Tin = /o(1500) and 



r,„^,, « 28 MeV, (66) 



r,„^,, « 26 MeV, (67) 



if (Till = /o(1710). The experiment gives for the decay of an into pions a 
value lying within the range 30 - 70 MeV . 



Now let us proceed with decays of am. The process din tt+tt is given 
by the amplitude 

that consists of two parts. The first part represents the contribution from 
the pure glueball. This contribution is small (since it is proportional to the 
pion mass squared), and the process is determined by the second part that 
describes the decay of the quarkonium component. As a result, the width of 
the decay am tttt if crni = /o(1500) is 

Tan,-^^^ = l^<rn,^n+n- « 14 McV, (69) 

and, if am = /o(1710), 

r.,„^.. « 8 MeV. (70) 

In the case of KK channels, the contribution of the pure glueball is also 
proportional to the kaon mass squared, and is rather large as compared to 
the pions case. The amplitude of the decay am K^K consists of three 
parts 
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where the pure glueball decay into K^K is given by the amphtude 

Xc 

The quarkonium contributions are 

\ 2 \FkI niu + ms J 

where and Fx are the pion and kaon weak decay constants, respectively, 
and Fg = nis/igsV^)- In the case when am is /o(1500), we have 

^aui^KK = ^aui^K+K- + ^aui^KORo = 2r^^^^^K+K- ~ 29 McV, (75) 

and in the other case {am = /o(1710)) 

^an^^KK « 60 McV. (76) 

The amphtude of the decay of am into 7777 and ijij' can also be consid- 
ered in the same manner. The only complication is the singlet-octet mixing 
in the pseudoscalar sector and additional vertices coming from ALan- The 
corresponding amplitude is 

^.ni^ryry = 2^^^6^v„i + 2^" siu^ Ob^^^^^^ + +2 A' cos^ ^6^^^^^ + 2Af sin^ eb^>,^^^, 

(77) 



(78) 
(79) 





M2 


A9 — 




rjrj 


Xc 


/tan _ 
^(f) — 






Xc 



where 6 = 6 — 60, with 6 being the singlet-octet mixing angle in the pseu- 
doscalar channel, 9 « —19° and 6*0 the ideal mixing angle, tan 6*0 = 1/a/2. 



The decay widths thereby are: if am = /o(1500), 
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r,„,_,, « 25 MeV, (80) 

and if am = /o(1710), 

r,„,_,, « 43 MeV. (81) 

For the decay of am into 7777', we have the following amplitude 

Aa,u^vv' = sin 2^6,„,„, + A' sin 2^6,^,„, - sin 2^6^v„:. (82) 

The direct decay of a bare glueball into 7777' is absent here. This process 
occurs only due to the mixing of the glueball and scalar isoscalar quarkonia 
and the anomaly contribution. The decay widths follows: 

Tain^W ~ 10 MeV, (83) 



for am = /o(1500), and 



Tain^W - 30 MeV. (84) 



for am = /o(1710). The estimate for the decay /o(1500) into 7/7/' is very rough, 
because the decay is allowed only due to a finite width of the resonance as its 
mass lies a little bit below the 7/7/' threshold. The calculation is made for the 
mass of /o(1500) plus its half-width. For /o(1710), we have a more reliable 
estimate, since its mass is large enough for the decay to be possible. 

Up to this moment we considered only decays into a pair of mesons. For 
the state crm, there is a possibility to decay into 4 pions. This decay can 
occur through intermediate a (/o(400 — 1200)) resonance. 

The decay through the cr-resonance can be represented as two processes: 
with two resonances crm aa ^ Att and one resonance crm a27T Att. 
The vertices determining these decays follow from Lagrangian (^B]). The 
decay of a glueball into two a is given by the amplitude: 

+2^-6^,.,„62^,^, (85) 
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where A^^^ is the pure glueball amphtudef^ 



M2 

« (86) 

Xc 

and the anomaly amphtude A^^ coming from ALan is 

AT ^ (87) 

The total amplitude describing the decay into 4 pions through two a- 
resonances is 

(A(^(S12)A^(S34) + 
+A(S14)A(S23)), (88) 

where the function Acr(s) appears due to the resonant structure of the pro- 
cesses 

A,(s) = {s- Ml + iM,^r,^y\ (89) 
where F^r is the decay width of the ai resonance (see ([63|)). This function 



depends on an invariant mass squared Sij defined as follows 

s,j = {k, + kjf, (z,j = l,...,4). (90) 
Here i and j enumerate the momenta ki of pions 7r+(/ci), 7r^(/c2), 7r+(A;3), and 

7r"(/C4). 

Now let us consider the decay into 47r through one cr-resonance. The 
amplitude describing this process is as follows: 

The glueball amplitude is 



•^^ To obtain an approximate estimate for the glueball contribution, we used the mass of state before 
diagonalization (see the term with cTg ^ in d52 
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and the quarkonium: 



K.. = (92) 

Xc 



^a:„^a2. = '^qIZ- (93) 



The gluebah contribution prevails over the quarkonium one in magnitude 
and is opposite in sign. 

The amphtude describing the decay am 27r~^27r~ through one a- 
resonance is: 

(A^(S12)+A^(S34) + 

+A,(si4)+A.(s23)), (94) 

The total amplitude of the decay into 27r+27r" via cr-resonances is obtained 
as a cumulative contribution from both one and two intermediate a mesons: 

^criii^27r+27r- = ^criii^crcr^27r+27r- + ^0-iii^cr27r^27r+27r- • (95) 

The amplitude describing the decay into 27r^7r+7r~ has the following form: 

^o'iii^27r°7r+7r~ ^(7iii^(7(7^27r°7r+7r~ ~l~ ^criii^cr27r^27r°7r+7r~ j (96) 

where 

Act(si2)A^(s34), (97) 

-2A^^^^^^2nAa^2A^a{si2) + A^(s34)). (98) 

In this case, ki and k2 are momenta of the two tt'^, and su is their invariant 
mass squared. Indices 3 and 4 stand for 7r+ and 7r~, respectively. The 
amplitude A(j^27r0 is equal to O.SAo-^tt+tt- 
In the case of the decay into Air^, we have 

^0-111^477° — ^iTiii^iTO-^47rO ~l~ ^criii^cr27r^47r0 5 (99) 

where 
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A^„i^^^_47rO = 4^ani^2a-4^^2^o(^a(si2)A^(s34) + A^(si3) A^(s24) + 

+ A,(S14)A,(S23)), (100) 
^criii^cr27r^27r"7r+7r- = " 274crjjj^cr27r^cr^27r0 ( Acr(si2) + Acr(si3) + Ao-(si4) + 

+ A^(S23) + A^(S24) + A^(S34)). (101) 

Let us give numerical estimates for these decay modes. The decay width 
of the glueball into four particles is calculated using the prescript given in 

M 

1 



r47r — 



64(27r)6rM2 



X 



'123 , 

X / ds 

^123 



123 



CTiii 



'ds 



'34 



ds:- 



12 / _ "<^34 

S34 Js. 



'23 



ds 



23 



23 



o'iii^47r I 



=, (102) 



-1 VA(S123,S12,M2)(1-C2: 

where A^r^^^^^j^ is the amplitude describing one of the processes discussed 
above, M^^^^ is the mass of am. The corresponding two-particle invariant 
masses are defined in (^) except for S123, the invariant mass of three pions 

S123 = {h + k2 + hf. (103) 

The cosine between the plane formed by 3-momenta ki, k2 and the plane 
formed by ka, k4 in the rest frame of three mesons (ki + k2 + ks = 0) is 
denoted by C. The limits of integration follows 



sr23 = 



5^23 = (^-m - M^)\ 



'12 



4^2, s^^{^i-M^] 



4 = + 



2si23 



[(S123 + Ml - Su){Ml^ - Ml - Si23)± 



±VA(S123, S12, M2)A(M2 5^23, M^) 



4 = 2Ml + -- [Si2(si23 - Ml - Su) ± 

2si2 



± VA(S12, Ml M2)A(S123, S12, M2) 



(104) 



Formula ( |iU2| ) is similar to that given in []3T|], however, we used here different 
kinematic variables. As a result, we obtain for the decay into 47r: 
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ra,u^2n+2n- « 2.2MeV, r,„^^2.o,+,- « 1.2MeV, r,„^^4.o « O.lMeV. 

(105) 

The total width is 

« 3.5MeV (106) 
and, in the other case {am = /o(1710)), 

r.,„^2.+2.- « 6MeV, r,„^^2.o,+,- « 3.3MeV r,„^^4.o « O.SMeV. (107) 
The total width is 

« lOMeV. (108) 

As one can see, these values are very small. 

The other possibility of the state am to decay into 4 pions is to produce 
two intermediate p-resonances {am ^ 2p ^ 47r). Contrary to the decay 
through scalar resonances, where strong compensations take place, in the 
process with p-resonances, no compensation occurs, and it turns out that 
the decay through p determines the most part of the decay width of crin- 

To calculate the amplitude describing the process am 2/?, we need a 
piece of the Lagrangian with /o-meson fields. Although we did not consider 
vector mesons in the source Lagrangian, an extended version of NJL model 
T3|] contains the vector and axial-vector fields. Taking the mass term for p 



mesons from |T3j and including dilaton fields into it according to the principle 



M2 /^^^2 
p 



of scale invariance, we obtain: 

(f]\2p;p-^+pX), (109) 

Vac/ 

where Mp = 770 MeV is the /o-meson mass. From this, we derive the vertex 
describing the decay am PP- 



M 



2 



^b^'ana'{2ptP~,+P>l)- (110) 
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The decay of a p-meson into pions is described by the following amplitude: 

9p{pi-nY- (111) 

where Qp = 6.14 is the /o-meson decay constant, pi and p2 are the momenta 
of 7r+ and tt^. Finally, we come to the following formula for the amplitude 
of the process am p^p^ 27r+27r~: 

^aiii^pV^27r+2^- = — ^-^^^( (Sl3 + S24 " Su " S23) Ap(si2) Ap(s34) + 

Xc \ 

+ (si3 + S24 - Su - S34) Ap(si4) Ap(s23) j , (112) 

The function Ap(s) is the following: 

A,(s) = (s - + zM,r,)-i. (113) 

Here Tp = 150 MeV is the decay width of the p resonance. The decay into 
27r*^7r+7r^ occurs through a pair of charged /o-resonances: p^ and p^. The 
amplitude of this process is the same as for the decay with intermediate p^. 
The decay into 47r'^ cannot go via /o-resonances. 

In an extended N JL model [|r3| , there are no vertices describing the decay 



of a quarkonium into p mesons. As a result, only the glueball part determines 
the decay of am into 4 pions through p resonances unlike the case with a 
resonances. This leads to a large decay rate through p mesons (contrary to 
decays through a). 

Now let us give the numerical estimates for the decay into 4 pions. In the 
case, where am = /o(1500) we have: 

^ani^pp^27:+27r- « 50, r^„j^^^^27r"^+7r- ~ 90 McV, (114) 

with the total width: 

Kt^An « 140MeV. (115) 
In the other case {am = /o(1710)), 

^aui^pp^27r+2n- ~ 350MeV, I" aui^pP^27rH+7r~ ~ 650 McV, (116) 
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C^4. « IGeV. (117) 

Now we can estimate the total width of the state am. If am is identified 
with /o(1500) we have 

r;;°* « 220 MeV, (118) 



which is in qualitative agreement with the experimental value 112 MeV [[29 
and, in the other case {am = /o(1710)), 

« 1.2 GeV, (119) 

which is in great contradiction with experimental data. In the last case 
(/o(1710)) p mesons can show up as on-mass-shell decay products at large 
probability. The decay width is estimated as ~1 GeV. The absence of this 
decay mode in experimental observations is a reason that /o(1710) is not a 
glueball. 

Our estimates for the decay widths of the scalar meson states di, dn, and 
am are collected in Table 3. 



7. CONCLUSION 

In the approach presented here, we assume that (with the exception of 
the dilaton potential and the 't Hooft interaction) scale invariance holds for 
the effective Lagrangian before and after SBCS in the chiral limit. On the 
other hand, we take into account effects of scale invariance breaking that 
come from three sources: the terms with current quark masses, the dilaton 
potential reproducing the scale anomaly of QCD, and term Lan induced by 
gluon anomalies (see (0) in the Introduction). 

The scale invariance breaking that is connected with the term Lan was not 
taken into account in our previous paper This led to a small quarkonia- 

glueball mixing proportional to current quark masses, disappearing in the 



''^Note that there was wrong sign at the term in formula (43) that describes the quarkonia-glueball mixing, 
which led to incorrect estimates for the decay widths of the scalar glueball. 
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chiral limit. If the term ALan taken into account in ( |25|) the quarkonia- 
glueball mixing becomes much greater and does not disappear in the chiral 
limit, being proportional to constituent quark masses (quark condensates). 
This corresponds to results obtained from QCD in [03| . This contribution to 
the quarkonia-glueball mixing turns out to have decisive effect on the strong 
decay widths of scalar mesons. 

For the scalar meson states /o(400 — 1200) and /o(980), we obtain good 
agreement with experimental data Their decay widths are determined 
by quarkonium parts of decay amplitudes. 

Strong decays of the scalar meson state am ("glueball") are considered 
for two different masses: 1500 MeV and 1710 MeV. In the tttt channel, the 
contribution from quarkonia prevails over that from the glueball and thereby 
determines the decay rate. In the case of KK, ijij, ijij'^ and tttt channels, 
there are noticeable compensations among decay amplitudes. 

A similar situation with compensations takes place in the decay into 47r 
with intermediate cr-mesons. Here we have a strong compensation among 
the glueball and quarkonia contributions. But there is a possibility for the 
state am to decay through /o-resonances. In this case, as the quarkonium 
component is absent, no compensation occurs, and this channel determines 
the most of the total decay width of am- 

We performed calculations for both candidates for the scalar glueball state: 
/o(1500) and /o(1710), and found that /o(1500) is rather the glueball. The 
main decay mode is that into 4 pions. The decay rate into a pair of kaons 
is next by order of magnitude and is followed by the 7777, 7777', and tttt decay 
modes. 

The total width of the third scalar isoscalar state is estimated to be about 
220 MeV for M^j„ = 1500 MeV and 1.2 GeV for M^j„ = 1710 MeV. The 
experimental width of /o(1500) is 112 MeV and that of /o(1710) is 130 MeV. 
Unfortunately, the detailed data on the branching ratios of /o(1500) and 
/o(1710) are not reliable and are controversial [j2^ . 
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Our calculations are rather qualitative. However, they allow us to conclude 
that /o(1500) is a scalar glueball state, whereas /o(1710) is a quarkonium, 
for the following reasons: 1) The total decay width of the glueball in our 
model better fits to its experimental value if /o(1500) is assumed to be the 
glueball, rather than /o(1710). 2)As it follows from our calculations, the 
main decay mode of the scalar glueball is that into four pions. This is true 
for the state /o(1500). A decay of /o(1710) into four pions, however, was not 
seen in experiment. 3) Moreover, a direct decay of into a pair of p mesons 
on their mass-shell is possible for a scalar glueball with the mass about 1.7 
GeV. It has also not been observed. Our conclusion concerning the nature 
of /o(1710) as a quarkonium state is in agreement with the conclusion made 
in our papers []T^ . 

We are going to use this approach in our future work for describing both 
glueballs and ground and radially excited scalar meson nonets which lie it 
the energy interval from 0.4 to 1.71 GeV. 
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8. TABLE CAPTIONS 



Table The masses of physical scalar meson states ai, an, am and the values of 
the parameters Xc, Xo, bag constant B, and (bare) glueball mass nig (in 
MeV) for two cases: 1) M^„j = 1500 MeV and 2) M^„j = 1710 MeV 

Table 0. Elements of the matrix b, describing mixing in the scalar isoscalar sector. 

The upper table refers to the case am = /o(1500), the lower one to the 
case am = /o(1710) 

Table The partial and total decay widths (in MeV) of the scalar meson states 
/o(400- 1200), /o(980) and of the gluebaU for two cases: am = /o(1500) 
and am = /o(1710), and experimental values of decay widths of /o(1500) 
and/o(1710) 



32 



9. TABLES 



TABLE 1. 









C^III 


Xc 


Xo 


B, GeV^ 


M, 


I 

II 


400 
400 


1100 
1100 


1500 
1710 


206 
180 


190 
166 


0.009 
0.009 


1447 
1665 


TABLE 2. 








o-i 






<7"II 




C^III 


V 

x' 






0.939 
-0.214 
-0.270 






0.240 
0.968 
0.067 




0.247 
-0.128 
0.960 










o-i 










<7III 


< 

x' 






0.948 
-0.216 
-0.233 






0.232 
0.971 
0.047 




0.216 
-0.099 
0.971 



TABLE 3. 





F 

■■- TTTT 


^KK 






r47r 


Ttot 


T-iexp 
tot 


/o(400 - 1200) 


820 










820 


600-1200 


/o(980) 


28 










28 


40-100 


/o(1500) 


14 


28 


25 


13 


~ 140 


~ 220 


112 


./o(-l-710) 


8 


GO 


13 


31 


- 1000 


- 1100 


130 
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